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Abstract. The use of a method of spectra disentangling for telluric lines is explained in detail, with a particular emphasis 
on high-precision radial-velocity measurements for the search for extrasolar planets. New improvements to the method are 
introduced. 
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1. Introduction 

Ground-based spectroscopy of astrophysical objects is influ- 
enced by telluric lines, which are often pronounced especially 
in the red and infrared region. These lines have to be taken into 
account for numerous reasons, the most important of which are: 
(i) the need to clean the stellar spectra of this additional absorp- 
tion, (ii) the possibility to check or to improve the definition of 
the wavelength scale in the spectra (see, e.g., Griffin and Griffin 
1973, Greimel and Yang 1999) and (iii) a study of the Earth's 
atmosphere (cf. Griffin 2005). While in the first case the telluric 
absorption is harmful, in the other two cases it is useful. For 
the purpose (ii) not only are the natural atmospheric lines used 
but artificially generated lines from hydrogen-fluoride, iodine 
or other absorption cells (cf. e.g. Campbell 1983, Butler et al. 
1996) can also be imprinted on the observed spectra. Because 
the problems of the data processing methods discussed here for 
absorption cells and atmospheric lines are in principle equiva- 
lent, we shall refer to both cases as the "telluric lines". 

In the above items it is necessary to accurately separate the 
telluric lines from the ('stellar') spectrum of the observed back- 
ground object. Several methods have been developed for this 
purpose with different advantages and limitations. In principle, 
the telluric lines can be distinguished from the stellar lines us- 
ing a proper model of their wavelengths, ratios of strengths and 
line profiles (cf. e.g. Adelman et al. 1996) or by their variabil- 
ity in wavelength (due to the annual motion of the Earth) or in 
strength (due to variations of the air mass and humidity). 

The method of disentangling spectra of multiple stars 
(Simon and Sturm 1994, Hadrava 1995, PI hereafter) gener- 
alized for variable line-strengths (Hadrava 1997, P2 hereafter) 
is suitable to separate telluric lines. This possibility has been 
incorporated in the code KOREL for Fourier disentangling 
(Hadrava 2004b, P3 hereafter) and has been widely used (e.g. 
Harmanec et al. 1997 etc.). 

Hrudkova and Harmanec (2005, HH hereafter) applied 
the KOREL disentangling of telluric lines to a Boo and ad- 
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vertised its usefulness for studies of late-type single stars. 
Unfortunately, they had not properly taken into account expla- 
nations of the method (Hadrava 2004a or P3) and give some 
misleading information. To prevent further confusion we give 
the explanation more explicitly here. In addition, some new 
possibilities and improvements of the method are presented. 

2. Disentangling of telluric lines using KOREL 

The disentangling of spectra is made possible by the varia- 
tions of mutual Doppler shifts, and thus works equally for dif- 
ferent components of a multiple stellar system or for telluric 
lines moving annually with respect to single or multiple stars. 
Telluric, interstellar and binary spectra can be disentangled si- 
multaneously. However, because the telluric lines also vary in 
strength, their disentangling became possible in practice only 
after a generalization enabling one to take into account that 
line-strength variation (cf. P2). On the other hand, this vari- 
ability in strength is a helpful feature in disentangling telluric 
lines. 

In the first applications, the orbital parameters of telluric 
lines were converged by the disentangling method. However, 
because these parameters are given by the annual motion of the 
Earth and by the ecliptic coordinates of the observed object, 
which are known, it is advantageous to use the constraints on 
these parameters as facilitated by the code PREKOR (cf. P3 for 
details). The observed stellar spectra are usually reduced to the 
heliocentric wavelength scale. In this case the fictitious centre 
of mass of the Earth's relative motion (i.e. orbit No. 3 if the 
telluric lines are identified with component No. 5 - cf. Fig. 1 in 
PI) with respect to the star or stellar system is put close to the 
star, imposing a very small if -velocity of the star, and a pro- 
portionally small value of the mass ratio q% = TOEarth/?ri s tar 
is chosen to give the proper amplitude i^Earth of the telluric 
lines. The value K stm — lm/s is preset and the correspond- 
ing value of <73 is calculated from the coordinates of the star by 
the code PREKOR, but users may rescale both these numbers 
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655.0 660.0 A[nm] 

Fig. 1. The standard KOREL output for three spectra of r Boo 
(upper three lines) taken on three consecutive nights disen- 
tangled as stellar (the fourth line) and telluric spectra (bottom 
line). 



from PREKOR to keep K stal negligible if a higher accuracy is 
desirable. This can be necessary especially in applications that 
search for extrasolar planets or other low-mass companions, 
when the lines of an absorption cell can also be superimposed 
on the stellar spectrum to provide a more precise wavelength 
scale (as mentioned in Section [Q. However, for such studies 
requiring high precision, the Keplerian orbit is an insufficient 
approximation for the velocity of the Earth and a more accurate 
approximation, which takes into account corrections for other 
planets in the Solar system, has to be involved. For this purpose 
the PREKOR code includes the subroutine kindly provided by 
C. Ron and J. Vondrak (Ron and Vondrak 1986) which calcu- 
lates the vector of the Earth's velocity with a precision of 0.02 
m/s. The disentangling can be then performed with the veloci- 
ties of the telluric iodine lines unbound to a Keplerian orbit (as 
described in P3) but fixed by this more precise model. Some 
differences between radial velocities of telluric lines in indi- 
vidual exposures found by KOREL and the Ron and Vondrak 
model can be caused by instrumental or data-processing er- 
rors in determining the wavelength scale. The whole spectrum 
can be then shifted by a new run of PREKOR. In principle, 
a version of the disentangling method with self-adjustment of 
wavelengths could be developed for cases with well-defined 
telluric/iodine lines in all spectra. 

An alternative possibility is to disentangle the spectra in in- 
strumental, instead of in heliocentric, wavelengths, to set the 
A-velocity of the whole observed stellar system equal to the 
corresponding projection of the Earth's velocity in the eclip- 
tic on the line of sight, and to formally set the mass ratio 
WEarth/wstar very large to decrease the A"- velocity of telluric 
lines safely below the resolution capability of the spectra. 



3. The role of line-strength variability 

It is claimed in HH that "the disadvantage of the method is that 
one needs a number of spectra of the star in question, secured at 
different months of the year" to disentangle telluric lines from 
spectra of a single star. It is not necessarily the case, as can 
be demonstrated by the example shown in Fig. ^ The three 
Ondfejov RETICON spectra in the region of H Q line of the star 
t Boo, known for spectroscopic observations of its planet, were 
taken within three days. To exclude the possibility that the dis- 
entangling of telluric lines is due either to the annual change 
of heliocentric correction of the radial velocity (which is about 
0.6 km/s between the first and the third exposure) or to the mo- 
tion of the star itself (which has a period of about 3 days and 
is detectable by disentangling a larger set of spectra), the times 
of all three exposures were set equal to the time of the middle 
one. Nonetheless, the disentangling is satisfactory. This is ow- 
ing to the differences in strengths of the telluric lines, which 
are obviously weakest in the second and strongest in the third 
spectrum. The difference between two chosen spectra with dif- 
ferent strengths determines the shape of the telluric spectrum 
Iteii (in some reference scale), which must be then subtracted 
from each exposure I i, s (t) with a proper multiplicative line- 
strength factor s(t) to get the stellar spectrum Iq, 

h = Iobs{t) - s(t)I teU . (1) 

The solution of factors s(t) is not unique, because ks(t) (k > 
0) gives the same solution with the rescaled telluric spectrum 
Iteii/k and s(t) — As fits the observed spectra equally by the 
stellar spectrum Iq + Aslteii- The strengths s(t) thus cannot 
be solved directly from equation (7) in P2 or Eq. (30) in P3, 
which is singular in this case, but can be converged by a sim- 
plex method which finds some of the existing solutions. By 
fixing some of the strengths we can choose the proper solution, 
or it can be more easily found by KOREL using even relatively 
small seasonal differences in heliocentric radial-velocity cor- 
rections (either for stars close to the poles of the ecliptic or 
for data covering a short part of the year). An important con- 
dition which is worth satisfying during the observation is thus 
obtaining spectra at different air masses. (Our early experience 
confirms the seemingly paradoxial fact that the disentangling of 
telluric spectra significantly improves if spectra that are highly 
blended with telluric lines are included in the solution.) 

Telluric lines consist of several subsystems whose strength 
may vary independently. Typically, lines of O2, N2 and Ar de- 
pend mainly on the air mass, while lines of H2O depend also on 
the instantaneous weather (cf. Adelman et al. 1996) and the O3 
lines should be checked for long-term variations (Griffin 2005). 
Because these lines have identical radial velocities, they can be 
disentangled only by their variable strengths. In KOREL, it is 
possible by the above explained method to assign the annual 
motion of the Earth to both orbits No. 2 and 3 (and a negligible 
if-velocity to orbit No. 1), so that all the components in No. 
3, 4 and 5 can correspond to telluric lines (cf. Fig. 1 in PI). 
However, to separate the subsystems, good initial estimates of 
their strengths are needed, or spectra of at least some of them 
should be constrained by an accurate model, as explained be- 
low. An alternative method, introduced by Lynch (2005, cf. 
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Lynch and Polidan 1997) to separate spectral components that 
change their line-strengths independently, could be used for ex- 
posures with negligible differences in radial velocities between 
the stellar and telluric lines. This also confirms that the limita- 
tion claimed by HH is incorrect. Lynch's method is based on 
the SVD/PCA technique and indicates how many independent 
components (characterized by their strengths) are present in the 
data, while this number must be fixed in the present versions of 
disentangling. Combination of both methods can thus be useful 
for the problems in question. 

4. Non-linearity of telluric lines 

HH mentioned as a "principal problem when KOREL is 
used to disentangle telluric lines" the question of nonlinearity 
(Hensberge 2001). HH then incompletely rephrased the origi- 
nal answer given by the present author and ignored the recent, 
more thorough explanations (Hadrava 2004a or P3). The 'prob- 
lem' is that the telluric spectrum is not an additive, but a mul- 
tiplicative component in the observed spectra, i.e. the observed 
spectrum is (cf. Eq. (31) in P3) 



Iob s {x,t) =exp(-T(x,t))I (x,t) , 



(2) 



where Iq is the spectrum of the studied object as seen outside 
the Earth's atmosphere and r is the optical depth along the line 
of sight in the atmosphere. This can be expressed as a sum of 
Iq and a negative telluric component: 



hen = (er T - l)/ - -tI 



(3) 



which also depends on Iq. Consequently, the depth of the tel- 
luric line in Fig. ^ which falls close to the centre of the stellar 
H Q line is underestimated in the decomposed spectrum relative 
to the depths of telluric lines falling on the continuum. If the 
disentangled spectra are to be used for spectroscopic studies 
of the atmosphere of the Earth, the telluric spectrum must be 
divided by the mean stellar spectrum. Such a renormalisation 
has also to be applied if we want to use the telluric spectrum 
disentangled from one star as a template for constrained disen- 
tangling of another star, as will be explained later. The only dis- 
crepancy arises if the telluric line moves (either due to its own 
seasonal shifts or due to the orbital motion of the star) in the 
part of stellar line-profile with a steep slope, or if the strength 
of the line varies, e.g. in an eclipse of the star. The error in the 
residuals (O-C) is proportional to the depth of the telluric line 
multiplied only by the difference in the local depth of the stellar 
line from its mean value, not its whole depth. Possible deflec- 
tion of the decomposed spectrum is generally even smaller. 

In the simple case of a one-component stellar system (ei- 
ther a single star or a multiple system dominated by the spec- 
trum of one component) plus a telluric spectrum the problem 
of nonlinearity of the telluric component can be avoided by a 
logarithmic transformation of the (rectified) input spectra: 



I(t,x) -> l'{t,x) = 1 + lnJ = 

= 1 + ln/ (£, x) — r(t, x) , 



(4) 



where Iq is the stellar spectrum and t is the monochromatic op- 
tical depth of the Earth's atmosphere (cf. Eq. (31) in P3). The 



corresponding backward transformation should be applied to 
the disentangled spectra. This would not be appropriate to dis- 
entangling standard binaries (for which the spectra are additive, 
Iq = h + h + ■■■)■ However, in searching for a faint compo- 
nent I2 in the presence of strong telluric lines in the spectrum 
I\ of the dominant component, it can be advantageous to use 
the linearization InJo = ln/i + (I2 + ■■■)/h and to disentan- 
gle the spectra on a logarithmic scale of intensities. Naturaly, 
the normalization of disentangled telluric spectrum described 
in the previous paragraph must not be performed if the spectra 
are transformed into a logaritmic scale of intensity. Instead, the 
inverse transforms of both the telluric and the stellar spectra 
have to be applied. 



5. Solution instabilities and use of constraints 

Despite the solutions described above, some problems may ap- 
pear or increase when telluric lines are included in the disentan- 
gling. One example is the problem of waves in the disentangled 
spectra Ij(x), due to instabilities of the lowest Fourier modes 
Ij(y), for which the corresponding equations (Eqs. (5) in PI, 
(6) in P2, (20) or (28) in P3) of general form 



E 



JV 



E Wl A? (y> *«> p) (y,ti, p) 



.1=1 



N 



^2wiI(y,ti)A* m (y,ti,p) 



(5) 



1=1 



are nearly singular (here / are the spectra observed at time ti, 
Aj are broadening functions for component j - cf. PI, P2, P3 
for details). These instabilities often increase with a new com- 
ponent such as a telluric spectrum, at least until the appropriate 
strengths or orbital parameters are well approximated. Possible 
errors in rectification or real non-Doppler variations of stellar 
spectra may force the solution to use its new freedom to fit these 
features instead of the true telluric spectrum. Two hints are rec- 
ommended in this case: (i) determine the strengths of telluric 
lines from a spectral region where they are dominant (e.g. the 
short-wavelength region shown in Fig.Q and fix them in other 
regions at least until the other parameters are not iterated fur- 
ther, or (ii) use the fact that the telluric lines are mostly narrow, 
i.e. contained in higher Fourier modes, and filter out the lowest 
modes at least until the correct solution is approached in higher 
modes. 

Quite often some problems caused by the rectification, line- 
profile variability or an unresolved component of the spectrum 
persist in centres of deep stellar lines and a small contribution 
to those lines then appears in the telluric spectrum, which is 
otherwise reasonable (with the correct lines on a flat contin- 
uum) in the far wings of the stellar lines and in the stellar con- 
tinuum. In such a case, use of a recently developed 'constrained 
disentangling' procedure can be helpful. This will be described 
elsewhere (Hadrava 2005). For our present purpose it is neces- 
sary to know that the observed spectra can be approximated by 
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the form 

m 

I{x,t) = ^I j {x)*A j {x,t,p) + 

n 

+ Jj(x) * Aj(x,t,p) , (6) 

j=m+l 

where the first sum includes the component spectra Ij which 
are unknown and have to be isolated by the disentangling and 
the second one includes those spectra which are known and 
constrained by some templates Jj (all component spectra being 
Doppler shifted or broadened by broadening functions Aj, the 
parameters p of which may be also solved by disentangling). 
The shape of the telluric spectrum can thus be constrained in 
this sense, either removing the residuals of stellar lines from the 
preliminarily disentangled spectra by their additional rectifica- 
tion, or directly by using a template for the telluric spectrum 
obtained from solutions for other stars (with the abovemen- 
tioned renormalization). With these constrained telluric spec- 
tra, the solution may be forced to divide the previous artifacts 
between the different disentangled stellar spectra (and to give 
better mean line-profiles), or to put the discrepancies into resid- 
ual spectra, from where their nature can be judged. 
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